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Abstract 



We discuss the relic density of the Ughtest of the supersymmetric particles 
(LSP) in view of new cosmological data, which favour the concept of an accel- 
5^ \ crating Universe with a non- vanishing cosmological constant. The new bound 

on the Cold Dark Matter density, r^cDM^o ~ 0.22, puts stringent constraints on 
supersymmetry preferring low supersymmetry breaking scales, in sharp contrast 
to electroweak precision measurements favouring large supersymmetry breaking 
scales. Supersymmetric predictions are in agreement with cosmological data and 
electroweak precision data in the window of the parameter space: tuq < 200 GeV, 
300 GeV < Ml/2 < 400 GeV, putting bounds on sparticle masses, which may be 
evaded if m-mp < nirj^ < 1-2 m^gp. 



Recent observations of type la supernovae (SNIa) put new constraints on the cosmo- 
logical parameters. The data favour an almost flat and accelerating Universe, where the 
acceleration mainly is driven by a non-vanishing cosmological constant. 

There is a growing consensus that the anisotropy of the Cosmic Background Radiation 
(CBR) offers the best way to determine the curvature of the Universe and hence the total 
matter-energy density VLq [0. The data are consistent with a flat Universe, since VLq = 
1.0±0.2 |l|,@], and the radiation content of the matter-energy density, that is contribution 
coming from CBR and/or ultra relativistic neutrinos, is very small. Therefore the present 
matter-energy density can be decomposed basically into two components: the matter 
density ^m and the vacuum energy Q/^: 

fio = ^M + ^A • (1) 

There is supporting evidence, coming from many independent astrophysical obser- 
vations, that the matter density weighs ^m = 0.4 ± 0.1. Recently two groups, the 
Supernova Cosmology Project and the High-z Supernova Search Team 0, using dif- 
ferent methods of analysis, each found evidence for accelerated expansion, driven by a 

vacuum energy contribution: 
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^A = g^M + 3 ± g . (2) 

So, for r^M = 0.4 ± 0.1 this relation implies that the vacuum energy is non-vanishing, 
Qa = 0.85 ±0.2, value which is compatible with a flat Universe, as the anisotropy of CBR 
measurements indicate. Taking into account the fact that the baryonic contribution to 
the matter density is small, Qb = 0.05 ± 0.005, the values for matter energy density Qu 
result to a Cold Dark Matter (CDM) density ^cdm — 0.35 ± 0.1, which combined with 
more recent measurements [||,0 of the scaled Hubble parameter ho = 0.65 ± 0.05, result 
to small CDM relic densities: 

^cdmC- 0.15 ±0.07. (3) 

Such stringent bounds for the CDM relic density do affect supersymmetric predictions 
and can severely lower the limits of the effective supersymmetry breaking scale and hence 
the masses of the supersymmetric particles, as first emphasized in Ref . . The CDM 
relic density with non-vanishing cosmological constant in the framework of the Minimal 



Supersymmetric Standard Model (MSSM) has been the subject of Ref. [0. Using the 
recent cosmological data, gauge fermions are predicted to be within LHC reach. 

It is perhaps worth pointing out that while electroweak (EW) precision data are 
in perfect agreement with Standard Model (SM) predictions, and hence in agreement 
too with supersymmetric models characterized by a large supersymmetry breaking scale 
^susY m, the data on ^cdm ^o^ push Msusy to the opposite direction favouring small 
values of Msusy- Therefore EW precision data may not reconcile with the assumption 
that the lightest supersymmetric particle (LSP or x), is a candidate for CDM. 

In i?-parity conserving supersymmetric theories the LSP is stable, and for most of 
supersymmetric models is the lightest neutralino, which is a good candidate for the CDM 



particle 0. Many authors |p!0|-p5[| have calculated the relic neutralino density. In the 
early works, only the most important neutralino annihilation channels were considered, 
but later works [^,|2T[] included all annihilation channels. Also more refined calculations 
of thermal averages of cross sections were employed, which took into account threshold 
effects and integration over Breit-Wigner poles PB|P7[ |. 

Our study in this letter is based on the constrained supergravity (SUGRA) scenario, 
assuming universal boundary conditions for the soft breaking parameters at the unifica- 



tion scale MgutQ- It is also assumed that the EW symmetry is radiatively broken [^ 



Therefore the arbitrary parameters are: mo, M1/2, ^0 aiid tan/5. The absolute value 
of /i is determined from the minimization conditions of the one-loop corrected effective 
potential. These also determine the Higgs mixing parameter mg. The sign of fi is unde- 
termined in this procedure and in our analysis both signs of /i are considered. Therefore 
in this scheme the fi value as well as m'^ = Bfi are not inputs. 

The Boltzmann transport equation for the neutralino number density is: 

^ = -3Hn - {avrei) {n'{T) - nl{T)) , (4) 

where n is the number of x's per unit volume, no is their density in thermal equilibrium 

and H is the expansion rate of the Universe. Using the fact that the total entropy 

S = hT^ R^ is conserved and defining the quantity q{x) = ^gLL with x = T/m-^, this 

^We allow however for small deviations from the gauge coupling unification scenario. In this case the 
value of the strong coupling constant at the unification scale, defined as the point where the couplings 
ai and 02 meet, is different from ai_2- 



differential equation can be cast into a form suitable for numerical manipulations [|TI[|: 

dQ ^^,,^ / 2 2^ 



where 



^^ A(x) (g^ - go^) , (5) 

m - (^G^) '^' ^^ (MT) + ^/.'(T)) (a..e.) . (6) 

The functions g{T) and h{T) appearing in the equation above, are the effective energy 
and entropy degrees of freedom respectively, and they determine the Universe energy 
and entropy density through the relations: 

P{T) = ^T' 9{T) , s{T) = ^ T^ HT) . (7) 

Solving Eq. (||) we obtain the q{xo), where Xq = Tq/ui^ corresponds to today's Universe 
temperature Tq ~ 2.7 °K. Using that q = y"LL , fi^ = -^ and pcrit = g^J' , we 
determine the present value of neutralino relic density Q^. h"^ from the equation: 

fl-hl = ^ T^ Gn m^ h{xo) q{xo) (lOOkmsec"^ Mpc"^)"^ (8) 

Before solving the Boltzmann equation we need to calculate the effective degrees of 
freedom functions g(T) and h(T), as well as the thermally averaged cross section {avrei), 
which enter into Eq. (|^). Regarding the calculation of the functions g{T) and h{T), 
the content of the particles in equilibrium is different depending on the temperature 



T. In our analyses we use the expressions for g{T), hiT) as given in Ref. ||T9l. Next 



we turn to the calculation of (aVrei)- At this point we follow Ref. pO| and express the 
non-relativistic cross sections for the various annihilation processes x X "^ ^ ^) i^ terms 
of helicity amplitudes. We follow the standard treatment and ignore contributions of all 
channels, which are forbidden at zero relative velocity Vrei of the two annihilating x's. 
This approximation is not expected to invalidate significantly our results. So we consider 
only the contributions of channels with non-supersymmetric particles in the final state: 

gg, /r, W-^W-, ZZ, ZH, Zh, ZA, W^H^, HH, hh, Hh, AA, HA, hA, H+R- . 

q and / denote quarks and leptons respectively, while if, /i, A and H^ denote the 
heavy, light, pseudoscalar and charged Higgses respectively. In our analysis we have not 
studied neutralino-stau coannihilation effects |l^,|2^,|2^, which if included can lower the 



values of the neutralino relic density. Although important, coannihilation processes are 
of relevance for values of the parameters for which m^ < rrif^ ^1-2 m^, that is near the 
edge where x &^d fa are almost degenerate in mass. 

It is well known that the non-relativistic expansion in the relative velocity Vrei breaks 
down near thresholds or poles of the cross sections, and in these cases, results based on 
this expansion are unreliable. We locate the points of the parameter space of the MSSM, 
which result to pole and/or thresholds of the cross section, using "near pole" and "near 



threshold" criteria. The comparison of our results with those of other studies [|17|, ^ , 
which treat the problem of poles and thresholds in a more accurate manner, shows that 
they are in striking agreement. This occurs, at least, in regions of the parameter space 
of MSSM where this comparison is feasible. 

Knowing {aVrei), from the procedure outlined previously, and by calculating the func- 
tions g{T), h{T), h'{T) we can have the prefactor A(x) appearing in Eq. (^. At high 
temperatures, or same large values of x = T/m^ above the freeze-out temperature, the 
function q{x) approaches its equilibrium value qo{x) (see Eq. ^. A convenient and ac- 
curate method for solving the Boltzmann equation is the WKB approximation employed 
in Ref. [l^. As far as the scanning of the parameter space of the MSSM is concerned, 
we exclude points that are theoretically forbidden, such as those leading to breaking of 
lepton and/or color number, or points for which Landau poles are developed and so on. 
We also exclude points for which the LSP is not a neutralino, as well as points of the 
parameter space for which violation of the experimental bounds on sparticle masses is 



encountered. We use the bounds which are listed in Ref. ||2^. From these bounds we 
have found that the chargino mass bound turns out to be the most stringent one. Details 
on our calculation will be published elsewhere pO|] . 



Before embarking to discuss our physics results we should stress that in our scheme we 

have not made any approximation concerning the masses or couplings of sectors which 

are rather involved such as neutralini for instance, which are crucial for our analysis. 

Therefore we do not only consider regions of the parameter space in which the LSP is 

purely B (bino) or purely a Higgsino, but also regions where in general the LSP happens 

to be an admixture of the four available degrees of freedomQ 

^ The case of a Higgsino- like LSP has been pursued in Refs. p3,p3], where the dominant radiative 
corrections to neutrahno masses are considered. Analogous corrections to couplings of Higgsino-like 



For large values of the LSP mass many channels are open, but for small values 
{m^ < 40 GeV) only channels with fermions, with the exclusion of top quark, in the 
final state are contributing. In these processes the exchanged particles can be either a 
Z-boson and a Higgs in the s-channel, as well as a sfermion / in the t-channel. Higgs 
exchanges are suppressed by their small couplings to light fermions, and sfermion ex- 
changes are suppressed when their masses are large. Then the only term surviving, even 
for large values of squark or slepton masses, is the Z-boson exchange. However in the 
parameter region where the LSP is a bino, this is not coupled to the Z-boson resulting 
to very small cross sections enhancing dramatically the LSP relic density. Therefore in 
considerations in which the LSP is a light bino^, large squark or slepton masses are in- 
evitably excluded, since they lead to large relic densities. If one relaxes this assumption 
and considers regions of the parameter space in which the LSP is light but is not purely 
bino, heavy squarks or sleptons may be allowed. We shall return to this later. Therefore 
the possibility for heavy g or / in the sparticle spectrum still exists, at the expense of 
having a light LSP and one of the chargino states. 

We have scanned the parameter space for values of ttiq, M1/2, Aq up to 1 TeV and 
tan/3 from around 1.8 to 35 for both positive and negative values of /i. The top quark 
mass is taken 175 GeV. In Figure ||, we display a representative output in the (Mi/2,'n^o) 
plane for fixed values of Aq and tan/5. Although in the displayed figure only values 
for /i > are presented, in our analysis both signs of /i have been considered. In the 
displayed figures Aq takes the values Aq = and tan/5 = 5, 20. The five different grey 
tone regions met as we move from bottom left to right up, correspond to regions in which 
Q^ hi takes values in the intervals 0.00 - 0.08, 0.08 - 0.22, 0.22 - 0.35, 0.35 - 0.60 and 
0.60 — 1.00 respectively^. 

In the blanc area covering the right up region, the relic density is found to be larger 
than unity. In the area to the left of the figure, the chargino mass bound is violated. 
Whenever a cross appears it designates that we are near either a pole or a threshold 
according to the criteria given previously. In these cases the approximations used are 



neutralinos to Z and Higgs bosons are important and can increase the relic density by a factor of 5, in 
regions of parameter space where LSP is a high purity Higgsino state |23] . 

•^ This happens when |/i| ^ Mw, with Mi smah « Mw- 

*These regions have been chosen in accord with the new bounds on H,^ Hq quoted in the introduction, 
and less stringent bounds on the same quantity which have been previously used in other works. 



untrustworthy and no safe conclusions can be drawn. For low values of M1/2 crosses 
correspond to mainly poles, which are either Z-boson or light Higgs, while for higher 
values, where LSP is heavier and hence new channels are open, these correspond to 
thresholds. For M1/2 ~ 110 GeV (^ 90 GeV for /i < 0) the lightest of the chargini 
has a mass close to its experimental lower limit. This bound is violated for all points 
^1/2 < 100 GeV. The dark area at the bottom part of the figure, which occurs for low 
values of mo, wiq < 100 GeV, is excluded since it mainly includes points for which the 
LSP is not a neutralino. In a lesser extend some of these correspond to points, which are 
theoretically excluded in the sense that either radiative breaking of the EW symmetry 
does not occur and/or other unwanted minima, breaking color or lepton number, are 
developed. 

For fixed M1/2 > 150 GeV the relic density f^x ^0 increases as tuq increases, due 
to the fact that cross sections involving sfermion exchanges decrease. Thus the area 
filled by the points corresponding to Qj^ /ig < 0.22 concentrates to the left bottom of the 
figure {rriQ < 200 GeV). Also for fixed rriQ < 250 GeV the number of points for which 
f2^ hi < 0.22 decreases as M1/2 increases since an increase in M1/2, enlarges squark and 
slepton masses as well yielding smaller cross sections. If M1/2 is further increased, for 
fixed niQ, the LSP will eventually cease to be a neutralinoQ. This picture changes for small 
values Ml/2 ~ 100 GeV, (see for instance the first of Figures P in which case small values 
of Qy. Hq < 0.22 can occur even for large values of the parameter ttiq > 300 GeV. This 
happens because the LSP is not a bino in this region and Z-boson exchanges contribute 
substantially to the annihilation of LSP's into light fermions, as explained previously. 
Actually in this region \fi\ is not much larger than the EW scale M^y and the LSP, 
although mostly bino, contains a sizeable portion of Hi Higgsino, leading to sizeable 
cross sections. 

Therefore the neutralino relic density does not exclude heavy squark or slepton states, 

as long as LSP is light, containing substantial mixture of HiggsinosQ. In all cases where 

this takes place, the LSP turns out to be light rriy. ^ 40 GeV or smaller and thus the only 

''The region for which the LSP is not a neutrahno depends rather strongly on the parameters ^0 a-nd 
tan/3. ActuaUy for large values of these parameters the light stau becomes lighter than the neutralinos. 



^ These corridors of low M1/2 and large toq values have been also presented in Ref. |21|. The figures 



of that reference displaying values of relic density in the (M^ /2, mo) plane are in remarkable agreement 
with our corresponding figures. 



open channels are those involving light fermions in the final state. Then the annihilation 
of LSP's into neutrinos for instance, a channel which is always open, proceeds via the 
exchange of a Z-boson which is non-vanishing and dominates the reaction when mo is 
sufficiently large, due to the heaviness of sfermions. This puts a lower bound on {aVrei) 
and hence an upper bound on Vt^ Hq, which can be lower than 0.22 consistent with the 
upper experimental limits quoted in the introduction. On these grounds one would 
expect that by increasing mo, while keeping M1/2 ~ 100 GeV fixed, the relic density 
stays below its upper experimental limit. However, with the exception of a few cases, 
this is not so in all cases studied. In some of the cases, beyond a certain point along 
the tuq axis, Q^ h^ exceeds the value 0.22 and starts increasing monotonically. This is 
due to the fact that the parameter |/i| gets large again making LSP moving towards 
the bino region. In this case cross sections are suppressed, since the coupling of LSP to 
Z-boson are negligible and sfermions are quite massive, and therefore the relic density 
is enhanced. As a further remark, we have to point out that the cosmologically allowed 
corridors of low M1/2 and high mo values, which we have just discussed, are almost ruled 
out in view of newer data on chargino masses [Q , which push up the lower bound for 
the soft gaugino mass M1/2 to about 130 GeV. 

It is seen from Figures |l] that as tan (3 increases the points for which the LSP is not a 
neutralino increase in number. As said before, this is due to the fact that by increasing 
tan/3 the stau sfermion f^ becomes lighter, since masses of the sfermions of the third 
generation have a rather strong dependence on tan/5. For exactly the same reasons this 
is also the case when we increase the soft parameter Aq. For lack of space we do not 
display this case [ pO[] . 

In Figures 0,0J5 we plot representative outputs of the relic density as function of one 
of the parameters mo, M1/2, ^0 and tan/5 keeping the other three fixed. In Figure ^ 
the relic density is plotted against mo] the lines shown correspond to different values for 
Ml/2- Notice that at the point where each line starts the mass of the LSP is equal to 
that of tr. From this figure it is obvious the tendency to get acceptable values for the 
relic density, as long as mo and M1/2 are kept light. Also obvious is the fact that larger 
values of mo are obtained for large values of tan (3. 

In Figure ^ the relic density is shown as function of the soft parameter M1/2. Low 



values of M1/2 < 120 GeV (110 GeV) for tan/3 = 5 (20) are not allowed because of the 
chargino mass bound. It is clearly seen that low values of M1/2 are preferred. The upper 
bounds on M1/2, which are compatible with the cosmo logical data, are higher for low mo 
and large tan (3. In this figure crosses denote location of poles or thresholds. 

In the first of Figures ^ we display the relic density as a function of Aq. We have 
chosen tan/3 = 20. Lowering the value of tan/3, keeping mo,Mi/2 fixed, the curves 
shown move upwards away from the shaded stripe, which is cosmologically accepted. In 
the second figure fi^ h^ is plotted as a function of tan /3. One immediately notices the 
tendency for the relic density to decrease as tan/3 gets large. 

Scattered plots of LSP relic density are shown in Figure ^. The sample consists of 
4000 random points that cover the most interesting part of the parameter space, which is 
within the limits : 1.8 < tan/3 < 40, M1/2 < 1 TeV, \Ao\ < 500 GeV and niQ < 500 GeVQ 
From the given sample only those points which lead to relics less than 1.5 are shown in 
the figure. The experimental bounds discussed before, restrict by about 40% the values 
of the allowed points. Low values of M1/2 less than about 100 GeV are experimentally 
forbidden by chargino searches. The points shown are striken by a cross (x) when 
mo < 100 GeV, by a plus (+) when 100 GeV < mo < 200 GeV and by a diamond (o) 
when mo exceeds 200 GeV. It is obvious the tendency to have M1/2 < 400 GeV in the 
cosmologically interesting domain which lies in the stripe between the two lines at 0.08 
and 0.22. Actually except for a few isolated cases, all allowed points are accumulated 
between M1/2 ~ 110 — 380 GeV. Note that the values of the parameter mo in the allowed 
area are restricted to mainly rriQ < 200 GeV (crosses or pluses). Only a few points with 
large values m-o > 200 GeV are consistent with the recent cosmological data. 

EW precision data are in perfect agreement with the SM and therefore with super- 
symmetric extensions of it which are characterised by a large Msusy- In fact overall fits 
to the EW precision data show a preference towards large values of Msusy, in which 
case better fits (lower x^) are obtained 0. In the constrained MSSM lower bounds on 
Ml/2 can be established from the experimental value of the effective weak mixing angle 
sin^ 6'efr, measured in SLC and LEP experiments. If the combined SLC and LEP data are 



^ Higher values for mo are of relevance only for M1/2 ~ 100 GeV. Such low values for M1/2 are 
almost ruled out by recent data. Also since flj^ /i§ does not depend strongly on Aq for A/1/2 > 100 GeV, 
as it can be realised from the first of the Figures 0, it suffices to focus on values |Ao| < 500 GeV. 



used, Ml/2 is restricted to be larger than about 300 GeV (see Dedes et. al. in Ref. 00). 
This lower bound on M1/2 can be further increased by about 200 GeV, if in addition 
gauge coupling unification is assumed^. 

Following the previous discussion, in Figure ^ we have shaded in grey the region which 
is consistent with EW precision data. This starts at about 300 GeV and progressively 
becomes darker as we move to larger M1/2 values, where the SM limit is attained and 
better agreement with experimental data is obtained. We observe that only a few points, 
the majority of them being very close to 400 GeV, are compatible with both astrophysical 
and EW precision measurements. The values of M1/2 and itlq for which theoretical 
predictions agree with both data, are constrained within the regions 300 GeV < A'/i/2 < 
400 GeV and itiq < 200 GeV. Performing a more refined scanning in the parameter 
space with values of M1/2, "^o in the aforementioned range and varying the other two 
in the range 1.8 < tan/3 < 40, l^ol < 1 TeV, we have found the following lower and 
upper bounds on the masses of the LSP, the lighter of charginos, staus, stops and the 
light scalar Higgs: 



"^LSP 


: 113 GeV 


- 149 GeV, 


m^ 


: 209 GeV 


- 278 GeV, 


rrifn 


: 118 GeV 


- 168 GeV, 


T^k 


: 390 GeV 


- 720 GeV, 


rrihr, 


: 76 GeV - 


- 122 GeV. 



These refer to the case /x > 0. Analogous bounds are found for yU < 0. With the 
exception of the Higgs lower bound, which is increased by about 20 GeV, the bounds 
for the remaining sparticles in the /x < case are almost unchanged. Since we have 
neglected coannihilation effects, the conclusions reached are actually valid outside the 
stripe m^ < nirj^ ^1-2 Tn^. Inside this band x~ ^ coannihilations dominate the cross 
sections, decreasing x I'elic densities leaving corridors of opportunity to high M1/2 and 
mo values as emphasized in other studies |2^ . 



^SLC data alone leave more freedom by allowing for lower M1/2 values. On the contrary LEP data 
alone prefer large values of M1/2. 

^However since this result depends sensitively on possible existence of High Energy thresholds we do 
not impose such a strong lower bound on Afi/2 values. 



In conclusion we can say that the new stringent bound on the matter rehc density 
Qy. /iq < 0.22 extracted from recent data, prefer low supersymmetry breaking scales and 
hence a light sparticle spectrum, provided rrif^ > 1.2 rriy.. This behaviour is exactly 
the opposite of what happens in the EW precision measurements physics. Reconciling 
LEP/SLC data with recent cosmological measurements restricts the parameters of the 
constrained minimal supersymmetry, yielding values of some of the sparticle masses not 
far away from their lower experimental limits. 
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A =0, tan/3=5, /j,>0 



A =0, tan/3=20, yLi>0 




Figure 1: The LSP relic density Q^^ h^ in the {mQ,Mi/2) plane for given values 
of Aq, tan/? and sign of /x. Grey tone regions, from darker to lighter, designate 
areas in which the LSP relic density takes values in the intervals: 0.00 — 0.08, 
0.08 - 0.22, 0.22 - 0.35, 0.35 - 0.60 and 0.60 - 1.00 respectively. In the 
blanc area ^^ ^o > -'-•O- "^^^ dark area corresponds to points for which the 
LSP is not a neutralino. The area which is excluded by chargino searches, 
ttIq > 66 GeV, is also shown. Crosses denote points for which thresholds or 
poles are encountered. 
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Figure 2: The relic density as function of ttiq for fixed values of Aq, tan/3 
shown in the figure. The solid line corresponds to M1/2 = 120 GeV (110 GeV) 
for tan/9 = 5 (20), the lowest allowed values by the chargino searches (see 
Figure |^). The dashed and dot-dashed lines correspond to M1/2 = 200 and 
400 GeV respectively. The grey area corresponds to Q^ h^ = 0.15 ± 0.07. 
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Figure 3: The relic density as function of M1/2 for fixed values of Aq, tan/3 
shown in the figure. The solid and dashed lines correspond to values rriQ = 
150 GeV and 200 GeV respectively. 
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Figure 4: The relic density as function of Aq and tan/5 for fixed values of the 
remaining parameters. The solid, dashed and dot-dashed lines correspond to 
Ml/2 = 120, 200 and 400 GeV respectively. 
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Figure 5: Scattered plot of the relic density versus M1/2 from a sample of 4000 
random points in the parameter space. Only the points with f2^ h^ less than 
1.5 are shown. The grey tone region within the cosmologically allowed stripe 
designates the region which agrees with EW precision data. 
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